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Abstract. It is shown that ambipolar diffusion as a toy 
nonlinearity leads to very similar behaviour of large scale 
turbulent dynamos as full MHD. This is demonstrated us- 
ing both direct simulations in a periodic box and a closure 
model for the magnetic correlation functions applicable to 
infinite space. Large scale fields develop via a nonlocal in- 
verse cascade as described by the alpha-effect. However, 
because magnetic helicity can only change on a resistive 
timescale, the time it takes to organize the field into large 
scales increases with magnetic Reynolds number. 



1. Ambipolar diffusion as a toy nonlinearity 

In this Letter we test and exploit the idea that the exact 
type of nonlinearity in the MHD equations is unessen- 
tial as far as the nature of large scale field generation is 
concerned. At first glance this may seem rather surpris- 
ing, especially if one pictures large scale field generation 
as the result of an inverse cascade process (Frisch et al. 
1975, Pouquet et al. 1976). Like the direct cascade in Kol- 
mogorov turbulence, the inverse cascade is accomplished 
by nonlinear interactions, suggesting that nonlinearity is 
important. However, a special type of inverse cascade is 
the strongly nonlocal inverse cascade process, which is 
usually referred to as the 'alpha-effect'; see Moffatt (1978) 
and Krause & Radler (1980). This effect exists already in 
linear (kinematic) theory. 

Until recently it was unclear which, if any, of the 
two effects (inverse cascade in the local sense or the a- 
effect) played the dominant role in large scale field gen- 
eration as seen in simulations (e.g. Glatzmaicr & Roberts 
1995, Brandenburg et al. 1995, Ziegler & Riidiger 2000) or 
in astrophysical bodies (stars, galaxies, accretion discs). 
A strong indication that it is actually the a-effect (i.e. 
the strongly nonlocal inverse cascade) that is responsi- 
ble for large scale field generation, comes from detailed 
analysis of recent three-dimensional simulations of forced 
isotropic non-mirror symmetric turbulence (Brandenburg 
2000, hereafter B2000). In those simulations a strong and 
nearly force-free magnetic field was produced, and most 



of the energy supply to this field was found to come from 
the forcing scale of the turbulence. 

In the absence of nonlinearity, however, the field seen 
in the simulations of B2000 became quickly swamped by 
magnetic fields at smaller scales. In that sense a purely 
kinematic large scale turbulent dynamo is impossible*. Any 
hope for analytic progress is therefore slim. However, the 
model of Subramanian (1997, 1999) is an exception. Sub- 
ramanian (1997; hereafter S97) extended the kinematic 
models of of Kazantsev (1968) and Vainshtein & Kitchati- 
nov (1986) by including ambipolar diffusion (in the strong 
coupling approximation) as a nonlinearity. Under the com- 
mon assumption that the velocity is delta-correlated in 
time, S97 derived a nonlinear equation for the evolution 
of the correlation functions of magnetic field and mag- 
netic helicity. Although the models of Kazantsev (1968) 
and Novikov et al. (1983) are usually known to describe 
small-scale field generation, Subramanian (1999; hereafter 
S99) found that in the presence of fluid helicity there is 
the possibility of tunnelling of bound-states correspond- 
ing to small scales to unbounded states corresponding to 
large scale fields, which are force-free. 

In this Letter we present numerical solutions to the 
closure model of S99. We stress that we do not advocate 
ambipolar diffusion (AD) as being dominant over the usual 
feedback from the Lorentz force in the momentum equa- 
tion. Instead, our motivation is to establish a useful toy 
model to study effects of nonlinearity in dynamos. Our 
numerical solutions may provide guidance for further an- 
alytic treatment of these equations in parameter regimes 
otherwise inaccessible. We begin however by considering 
first solutions of the fully three-dimensional MHD equa- 
tions in a periodic box using AD as the only nonlinearity. 



2. Box simulations for a finite system 

In this section we adopt the MHD equations for an isother- 
mal compressible gas, driven by a given body force /, in 
the presence of AD, but ignoring the Lorentz force 
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where D/Di = d/dt+u-'V is the advective derivative, B — 
V x A is the magnetic field, J = V x J3//io is the current 
density, and / is the random forcing function as specified 
in B2000. The nonlinear drift velocity uu due to AD can 
be written as «d = aJxB. We use nondimensional units 
where c s = k\ = po = Mo = 1- Here, c s is the sound speed, 
k\ the smallest wavenumber of the box (so its size is 27r), 
Po is the mean density, and /j,q is the vacuum permeability. 
Since AD is the only nonlinearity in Eq. (||) we can always 
normalize B such that a = 1. 

The model presented here is similar to Run 3 of B2000, 
where /i = rj = 2 x 1CP 3 . With a root-mean-square veloc- 
ity of around 0.3 the magnetic Reynolds number based on 
the size of the box is around 1000. The forcing wavenum- 
ber kf is chosen to be 5. In Fig. [I] we show a grey scale 
representation of a slice of the magnetic field and the cur- 
rent density at t = 337. Note the presence of a large scale 
magnetic field that varies in the z-direction. In Fig. || we 
show the spectra of magnetic and kinetic energies. The 
peak of magnetic energy at k = 1 shows the development 
of large scale magnetic fields. Further, the current density 
is concentrated into narrow filamentary structures, typical 
of AD (see Brandenburg & Zweibel 1994). 
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Fig. 1. Images of B y and J y in an arbitrarily chosen yz 
plane. 120 3 meshpoints, t — 337. 

Unfortunately, the severity of the diffusive timestep 
limit, St < 0.16fe 2 /?7AD, where ?7ad = aB 2 , prevented us 
from running much longer at high resolution (120 3 mesh- 
points). For 60 3 meshpoints this limit is unimportant, and 
so we were able to run until t = 900, a time when the large 
scale field was much more clearly defined. In the inset of 
Fig. 2, we show the evolution for such a case, but with a 
forcing at k{ = 10 (giving larger scale separation). Note 
again the peak of Em at k = 1 and also the suppression 
of magnetic field at the next smaller scale, corresponding 
to k > 2. Both these features are very similar to the mag- 
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Fig. 2. Spectra of magnetic energy (solid lines), kinetic 
energy (dash-dotted line), magnetic helicity (normalized 
by k/2; dashed line) for the run shown in Fig. [l| The 
inset shows spectra of a run with forcing at k = 10 and 
60 3 meshpoints for different times till t = 900. 

nctic field evolution in the case with full Lorentz force and 
without AD (Figs 3 and 17 of B2000). 

Our main conclusion from these results is first of all 
that large scale field generation works in spite of AD, con- 
trary to earlier suggestions that AD might suppress the 
large scale dynamo process (Kulsrud & Anderson 1992). 
Secondly, AD provides a nonlinear saturation mechanism 
for the magnetic field at all scales, except for the scale 
of the box, where a force-free field develops for which ud 
vanishes. Like in the simulations of B2000 this provides a 
'self-cleaning' mechanism, without which the field would 
be dominated by contributions from small scales. 

Having established the close similarity between models 
with AD versus full Lorentz force as nonlinearity, we now 
move on to discuss the nonlinear closure model of S99 with 
AD as a 'toy' nonlinearity. 

3. Closure model for an infinite system 

Under the assumptions that the velocity is delta- 
correlated in time and the magnetic field is a gaussian 
random field S97 derived equations for the longitudinal 
correlation function M(r, t) and the correlation function 
for magnetic helicity density, N(r,t). The velocity is rep- 
resented by a longitudinal correlation function T(r) and a 
correlation function for the kinetic helicity density, C(r). 
We change somewhat the notation of S99 and define the 
operators 
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so the closure equations can be written as 
M = 2D(t 1t DM) + 2GM + 4aH, 
N = -2i] T H + aM, 
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where H = —DDN is the correlation function of the cur- 
rent helicity, G = —DDT is the effective induction, 



'/T 



a (r) + 4aH(0,t) 
= r] + m (r) + 2aM(0,t) 
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are functions resembling the usual a-effect and the total 
magnetic diffusivity. Here cto(r) = — 2[C(0) — C(r)) and 
770 (r) = T(0) - T(r). Note that at large scales 



a(r 



\t(u-u) + Wah(J-B)/po (9) 



Vac = Vr(r -> 00) = \t(u 2 ) + ir A D<-B 2 )/MoPo, (10) 

where tad = 2apo- Expression (j^) is very similar to the a- 
suppression formula first found by Pouquet et al. (1976). 
Here a and tjt are scale dependent (i.e. they are largest 
on large scales) and, in addition, both are affected by AD. 

We construct T(r) and C(r) from an analytic approx- 
imation of the kinetic energy and helicity spectra, E^(k) 
and Hjt(k), respectively. Zero velocity at large scales 
means that Ex(k) ~ fc 4 for k — > 0. At some wavenum- 
ber fc = fcf the spectrum turns to a fc -5 / 3 Kolmogorov 
spectrum, followed by an exponential cutoff, so we take 



E K (k) 
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cxp(-fc/fc d ). 



(11) 



(fc/fcf) 17 /3 

We use parameters representative of the simulations of 
B2000, so E = 0.01, fcf = 5 and fc d = 25. Like in B2000 
we assume the turbulence fully helical, so Hk — 2fci?K 
(e.g. Moffatt 1978). The correlation functions T(r) and 
C(r) are then obtained via 



T(r) 



E(k)^^-dk = X(E(k)), 
kr 



(12) 



and C(r) = T(F(k))/4, where ji(x) = (sinx — x cos x)/x 2 
and r is the correlation time. (We use r = 4, representa- 
tive of the kinematic stage of Run 3 of B2000.) 

We solve Eqs. (||) and @ using second order finite 
differences and a third order time step on a uniform mesh 
in < x < L with up to 10,000 meshpoints and L = IOtt, 
which is large enough so that the outer boundary does 
not matter. In the absence of helicity, C = 0, and without 
nonlinearity, a = 0, we recover the model of Novikov et al. 
(1983). The critical magnetic Reynolds number based on 
the forcing scale is around 60. In the presence of helicity 
this critical Reynolds number decreases, confirming the 
general result that helicity promotes dynamo action (cf. 
Kim & Hughes 1997, S99). In the presence of nonlinearity 
the exponential growth of the magnetic field terminates 
when the magnetic energy becomes large. After that point 
the magnetic energy continues however to increase nearly 
linearly. Unlike the case of the periodic box (Sect. ||) the 
magnetic field can here extend to larger and larger scales; 
see Fig. |[ The corresponding magnetic energy spectra, 

f-L 
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Fig. 3. Evolution of magnetic correlation functions for 
different times. The correlation function of the magnetic 
helicity is shown in the inset. ?/ = 10~ 3 . 



10 







10 
£ 10 
10 
10 



10 

-1 



-2 
-3 
-4 



t= 100 



10 



10 



0.1 



1.0 
k 



10.0 



Fig. 4. Evolution of magnetic energy spectra. Note the 
propagation of magnetic helicity and energy to progres- 
sively larger scales. The fc -2 slope is given for orientation. 

are shown in Fig. |]. 

The resulting magnetic field is strongly helical and 
the magnetic helicity spectra (not shown) satisfy Hm ^ 
(2/k)Eu- The development of a helicity wave travelling 
towards smaller and smaller fc, as seen in Fig. is in agree- 
ment with the closure model of Pouquet et al. (1976). In 
the following we shall address the question of whether or 
not the growth of this large scale field depends on the mag- 
netic Reynolds number (as in B2000). We have checked 
that to a very good approximation the wavenumber of the 
peak is given by 



fcpeak(^) 



S)hoo(t). 



(14) 



This result is familiar from mean-field dynamo theory (see 
also S99) and is consistent with simulations (B2000, sec- 
tion 3.5). Note that here fc poa k decreases with time because 
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tends to a finite limit and r/oc increases. (This is not 
the case in the box calculations where fc pea k > 2-k/L.) 

4. Resistively limited growth on large scales 

In an unbounded system the magnetic helicity, {A ■ B) — 
6./V(0, t), can only change if there is microscopic magnetic 
diffusion and finite current helicity (J ■ B) = 6iJ(0,i), 



d{A- B)/dt = -2tj{J ■ B). 



(15) 



The closure model of S97 and S99 also satisfies this con- 
straint. (Note that ambipolar and/or turbulent diffusion 
do not enter!) As explained in B2000, this constraint limits 
the speed at which the large scale field can grow, but not 
its final amplitude. One way to relax this constraint is if 
there is a flux of helicity through open boundaries (Black- 
man & Field 2000, Kleeorin et al. 2000), which may be 
important in astrophysical bodies with boundaries. Here, 
however, we consider an infinite system. 

In Fig. H we show that, after some time t = i s , (J ■ B) 
reaches a finite value. This value increases somewhat as 
7/ is decreased. In all cases, however, (J ■ B) stays below 
(u: ■ u)(2t/ci), so that {a^l remains finite; see (Bh. A con- 
stant (J ■ B) implies that (A ■ B) grows linearly at a rate 
proportional to r\. However, since the large scale field is 
helical, and since most of the magnetic energy is by now 
(after t = t s ) in the large scales, the magnetic energy is 
proportional to (B 2 ) « fc poa k(A • B), and can therefore 
only continue to grow at a resistively limited rate, see 
Fig. |. 

5. Conclusions 

Our results have shown that ambipolar diffusion (AD) pro- 
vides a useful model for nonlinearity, enabling analytic (or 
semi-analytic) progress to be made in understanding non- 
linear dynamos. There are two key features that are shared 
both by this model and by the full MHD equations: (i) 
large scale fields are the result of a nonlocal inverse cas- 
cade as described by the a-effect, and (ii) after some initial 
saturation phase the large scale field continues to grow at 
a rate limited by magnetic diffusion. We reiterate that in 
astrophysical bodies the presence of open boundaries may 
relax the helicity constraint. Furthermore, the presence of 
large scale shear or differential rotation provides a means 
of amplifying toroidal magnetic fields quite independently 
of magnetic helicity, but this still requires poloidal fields 
for which the above conclusions hold. 
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Fig. 5. (a) Evolution of (J ■ B) for different values of ij. 
The corresponding value of is shown on the right hand 
side of the plot, (b) The evolution of magnetic energy for 
the same values of rj. 
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